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Abstract 

In a computer experiment the choice of suitable estimators to measure a physical quantity plays 
an important role. We propose a new direct route to determine estimators for observables which 
do not commute with the Hamiltonian. Our new route makes use of the Hellmann and Feynman 
theorem and in a diffusion Monte Carlo simulation it introduces a new bias to the measure due to 
the choice of the auxiliary function. This bias is independent from the usual one due to the choice 
of the trial wave function. We used our route to measure the radial distribution function of a spin 
one half Fermion fluid. 
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An important component of a computer experiment of a many particles system, a fluid, is the 
determination of suitable estimators to measure, through a statistical average, a given physical 
quantity, an observable. Whereas the average from different estimators must give the same result, 
the variance, the square of the statistical error, can be different for different estimators. We will 
denote with (0)f the measure of the physical observable O and with (. . .)y the statistical average 
over the probability distribution /. In this communication we use the word estimator to indicate 
the function O itself, unlike the more common use of the word to indicate the usual Monte Carlo 
estimator YJi=\ Oj/N of the average, where {(9, } is the set obtained evaluating O over a finite number 
N of points distributed according to /. This aspect of finding out different ways of calculating 
quantum properties in some ways resembles experimental physics. The theoretical concept may 
be perfectly well defined but it is up to the ingenuity of the experimentalist to find the best way of 
doing the measurement. Even what is meant by "best" is subject to debate. 

In ground state Monte Carlo simulations [HI Eh unlike classical Monte Carlo simulations [HI 
0, H] and path integral Monte Carlo simulations [g], one has to resort to the use of a trial wave 
function yj], While this is not a source of error, bias, in a diffusion Monte Carlo simulation 
|0] of a system of Bosons, it is for a system of Fermions, due to the sign problem [@]. Since this 
is always present in a Monte Carlo simulation of Fermions we will not consider any further when 
talking about the bias. 
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Another source of bias inevitably present in all three experiments, which we will not take into 
consideration in the following, is the finite size error. In the rest of the paper we will generally 
refer to the bias to indicate the error (neglecting the finite size error and the sign problem) that we 
make when defining different estimators of the same quantity not giving the same average. 

In a ground state Monte Carlo simulation, the energy has the zero- variance principle Q8|] : as 
the trial wave function approaches the exact ground state, the statistical error vanishes. In a dif- 
fusion Monte Carlo simulation of a system of Bosons the local energy of the trial wave function, 
E L (R) = /^(R), where R denotes a configuration of the system of particles and H is 

the Hamiltonian assumed to be real, is an unbiased estimator for the ground state. For Fermions, 
the ground state energy measurement is biased by the sign problem. For observables O which 
do not commute with the Hamiltonian, the local estimator, L (K) = [O v F(R)]/ v P(R), is inevitably 
biased by the choice of the trial wave function. A way to remedy to this bias can be the use of the 
forward walking method [§, 10] or the reptation quantum Monte Carlo method [11] to reach pure 
estimates. Otherwise this bias can be made of leading order S 2 , with 6 = <po - *P where (po is the 
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ground state wave function, introducing the extrapolated measure, O = 2(0 l)/ - (Ol)/^ where 
the first statistical average, the mixed measure, is over the diffusion Monte Carlo (DMC) station- 
ary probability distribution / and the second, the variational measure, over the variational Monte 
Carlo (VMC) probability distribution f vmc which can also be obtained as the stationary probability 
distribution of a DMC without branching IU2I1 . 

One may follow different routes to determine estimators such as the direct microscopic route, 
the virial route through the use of the virial theorem, or the thermodynamic route through the use 
of thermodynamic identities. In an unbiased experiment the different routes to the same observable 
must give the same average. 

In this communication we propose to use the Hellmann and Feynman theorem as a direct route 
for the determination of estimators in a diffusion Monte Carlo simulation. Some attempts in this 



direction have been tried before [|13l 11411 . The novelty of our approach, respect to Ref. I113ll . 



is a different definition of the correction to the variational measure, necessary in the diffusion 
experiment, and, respect to Ref. 014], the fact that the bias stemming from the sign problem does 
not exhaust all the bias due to the choice of the trial wave function. 

We start with the eigenvalue expression (H A - E A ) X ¥ A = for the ground state of the perturbed 
Hamiltonian H A = H + AO, take the derivative with respect to the parameter A, multiply on the 
right by the ground state at A = 0, (f> , and integrate over the particles configuration to get 



dR(p (H A -E A ) 



dE A dW 



dA dA 



W A 



Then we note that due to the Hermiticity of the Hamiltonian the left hand side vanishes at A 
so that we get further 



f dR<f> ¥ A 



dE A 



dA 



(1) 



A=0 



This relation holds only in the A —> limit unlike the more common form ||15|| which holds for 
any A. Given E A = f dR(p H Ax ¥ A / f dRcpo^ the "Hellmann and Feynman" (HF) measure in a 
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diffusion Monte Carlo experiment is then defined as follows 



— hf dE y 



A=0 



dA 

The a correction is 

#¥'(R) 



<0 L (R)>, + (AO a L (R)) f + (A(f L (R)) f . (2) 



AO l UR) = 



¥'(R) 



E L (R) 



»F'(R) 

— — . (3) 
¥(R) 



In a variational Monte Carlo experiment this term, usually, does not contribute to the average (with 
respect to f vmc oc *F 2 ) due to the Hermiticity of the Hamiltonian. We will then define a Hellmann 
and Feynman variational (HFv) estimator as HFv = Ol + AO". The correction is 

A(f L (R) = [E L (R) - E ] |^ , (4) 

where E = E A=0 is the ground state energy. It should be noticed that our correction differs by a 
factor 1 /2 from the zero-bias correction defined in Ref. [|13ll because these authors chose E A = 



J dWV A H Ax V A / JjR( v F 1 ) 2 right from the start. This correction is necessary in a diffusion Monte 
Carlo experiment not to bias the measure. The extrapolated Hellmann and Feynman measure will 
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then be O = 20 - (0 HFv )/ vmc . Both corrections a and to the local estimator depend on the 
auxiliary function, *P' = d x V A /dA\ / i = Q. Of course if, on the left hand side of Eq. ©, we had chosen 
v P /i=0 as the exact ground state wave function, <p , instead of the trial wave function, W, then both 
corrections would have vanished. When the trial wave function is sufficiently close to the exact 
ground state function a good approximation to the auxiliary function can be obtained from first 
order perturbation theory for A «: 1 . So the Hellmann and Feynman measure is affected by the 
new source of bias due to the choice of the auxiliary function which is independent from the bias 
due to the choice of the trial wave function. 

We applied the Hellmann and Feynman route to the measurement of the radial distribution 
function (RDF) of the Fermion fluid studied by Paziani Rial . This is a fluid of spin one-half 
particles interacting with a bare pair-potential v^(r) = erf (fir)/r immersed in a "neutralizing" 
background. The pair-potential depends on the parameter fx in such way that in the limit \x — » 
one recovers the ideal Fermi gas and in the limit // — » oo one finds the Jellium model. We chose 
this model because it allows to move continuously from a situation where the trial wave function 
coincides with the exact ground state, in the /i — > limit, to a situation where the correlations due 
to the particles interaction become important, in the opposite fi — » oo limit. 



We chose as auxiliary function W = QW, the first one of Toulouse et al. [17] (their Eq. (30)), 



r 2 K CdQ r 1 

Q^Ar, R) = -r-^ V <W,<W 7 -r-, ; , (5) 

%nVn (T n (T/ ^ ' J 4n r - r ; 



here <r and <x' denote the spin species, r = |r| the separation between two particles, r (J the separa- 
tion between particle i and j, (Tj the spin species of particle i, and dQ, r is the solid angle element of 
integration. The particles are in a recipient of volume V at a density n = n + + n_ = 1 /[47r(aor s ) 3 /3] 
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with do the Bohr radius, a = aor s the lengths unit, and the density of the spin <x parti- 
cles. With this choice the a correction partially cancels the histogram estimator I a>(r >(r, R) = 
Yii,j*id<T,o-i8<T',<Tj f S(r - rij)dQ.r/(4nVn a -n a ->), and one is left with a HFv estimator which goes to 
zero at large r. This is because the quantity (Al^^r, R))-^ = - f dyN ^(fVtVnQo-^ir, R) • dS/r 2 
equals minus one for all r with r e V, instead of zero as normally expected. This is ultimately re- 
lated to the behavior of the auxiliary function on the border of V N . The measure of the f3 correction 
also goes to zero at large r because one is left with a statistical average of a quantity proportional 
to E L (R) - Eq. The Hellmann and Feynman measure needs then to be shifted by +1. 

Our variational Monte Carlo experiments showed that in the variational measure the average 
of the histogram estimator agrees with the average of the HFv estimator within the square root of 
the variance of the average cr av = ^cr 2( KIN (here <t 2 is the variance, 7C the correlation time of 
the random walk, and N the number of Monte Carlo steps) and the two <x av are comparable. This 
is expected since the HFv estimator is defined exactly as in Ref. [[Till which correctly takes into 
account the definition of the HF estimator within a variational Monte Carlo simulation. In the fixed 
nodes diffusion experiment, where one has to add the f3 correction not to bias the average (note 
once again that this is defined by us as one half the zero-bias correction of Ref. II 1310 . the Hellmann 
and Feynman measure has an average in agreement with the one of the histogram estimator but 
the <r av increases. This is to be expected from the extensive nature of the f3 correction in which the 
energy appears. Of course the averages from the extrapolated Hellmann and Feynman measure 
and the extrapolated measure for the histogram estimator also agree. 

In the simulation for the Coulomb case, fi — > oo, we made extrapolations in time step and 
number of walkers for each value of r s . Given a relative precision 6 eo = Aeo/e x p , where eo = 
(E L )f/N, Ae is the statistical error on e > and e x p is the exchange energy, we set as our target 
relative precision 6 eQ = 10~ 2 %. The extrapolated values of the time step and number of walkers 
were then used for all other values of /i. We chose the trial wave function of the Bijl-Dingle-Jastrow 
11181 Il9l 12011 form as a prod uct of Slater determinants and a Jastrow factor. The pseudo-potential 
was chosen as in Ref. 112 ill . 3i, which is expected to give better results for Jellium. Comparison 
with the simulation of the unpolarized fluid at r s = 1 and fi = 1 with the pseudo potential of Ref. 
[|22ll . Sf\i f° r which the trial wave function becomes the exact ground state wave function in the 
fi — » limit, shows that the two extrapolated measures of the unlike histogram estimator differ one 
from the other by less than 7 x 10 -3 , the largest difference being at contact (see the inset of Fig. 
1). The use of more sophisticated trial wave functions, taking into account the effect of backflow 
and three-body correlations, is found to affect the measure by even less. In Tabled] we compare 
the contact values of the unlike RDF of the unpolarized fluid at various r s and /i from the measures 
of the histogram estimator and the HF measures. We see that there is disagreement between the 
measure from the histogram estimator and the HF measure only in the Coulomb n — > 00 case at 
r s = l,2. 

In conclusions we defined a Hellmann and Feynman estimator to measure a given physical 
property either in a variational Monte Carlo experiment and in a diffusion Monte Carlo experiment. 
Our definition coincides with the one of Ref. IU3I1 in the variational case but is different in the 
diffusion case. We proof tested our definitions on the calculation of the radial distribution function 
of a particular Fermion fluid. Our simulations showed that the bias is correctly accounted for 
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Table 1 : Contact values for the unlike RDF of the unpolarized fluid of Paziani Il6ll at various r s and /i from the mixed 
measure of the histogram estimator (hist) and the HF measure (HF) with the auxiliary function chosen as in Eq. Q, 
also reported are the two extrapolated measures (ext and HF-ext). The trial wave function used was of the Slater- 
Jastrow type with the Jastrow of Ref. lElll . J/2- The last column gives the error on the HF measure. 162 particles were 
used with 10 5 Monte Carlo steps. 
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in both kind of experiments but the variance increases in the diffusion experiment relative to the 
one of the histogram estimator. We believe it is still an open problem the one of determining the 
relationship between the choice of the auxiliary function and the variance of the Hellmann and 
Feynman measure. 

The idea for the work came from discussions with Saverio Moroni. I would also like to ac- 
knowledge the hospitality of the National Institute for Theoretical Physics (NITheP) of South 
Africa where the work was done. The simulations were carried out at the Center for High Per- 
formance Computing (CHPC), CSIR Campus, 15 Lower Hope St., Rosebank, Cape Town, South 
Africa. 
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Figure 1: Unlike RDF for the unpolarized fluid of Paziani IU6I1 at r s = 1 and fi — 1 with 162 particles. On the left 
panel the calculation with the Jastrow STi with various measures: variational histogram (variational) and variational 
HFv (HFv), mixed histogram (mixed) and HF (HF), and extrapolated histogram (extrapolated). On the right panel the 
calculation with the Jastrow J~2 with the histogram variational (variational JTzK mixed (mixed JTz), and extrapolated 
(extrapolated J2) measures. Also the extrapolated measure with the Jastrow J2 is compared with the extrapolated 
measure with the Jastrow J7i ■ In the inset is shown the difference between the histogram extrapolated measure of the 
calculation with Ji and the histogram extrapolated measure of the calculation with J'j. 10 5 Monte Carlo steps were 
used in the simulations. 
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